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Abstract 



For n G [—2,2] the 0{n) model on a random lattice has critical points to which a 



scaling behaviour characteristic of 2D gravity interacting with conformal matter fields 
with c G [— oo, 1] can be associated. Previously we have written down an exact solution 
of this model valid at any point in the coupling constant space and for any n. The 
solution was parametrized in terms of an auxiliary function. Here we determine the 
auxiliary function explicitly as a combination of ^-functions, thereby completing the 
solution of the model. Using our solution we investigate, for the simplest version of the 
model, hitherto unexplored regions of the parameter space. For example we determine 
in a closed form the eigenvalue density without any assumption of being close to or at 
a critical point. This gives a generalization of the Wigner semi-circle law to n 7^ 0. We 
also study the model for \n\ > 2. Both for n < —2 and n > 2 we find that the model 
is well defined in a certain region of the coupling constant space. For n < —2 we find 
no new critical points while for n > 2 we find new critical points at which the string 
susceptibility exponent •jstr takes the value +|. 
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1 Introduction 



For n G [—2,2] the 0{n) model on a random lattice has critical points to which 
a scaling behaviour characteristic of 2D gravity interacting with conformal matter 
fields with central charge c G [— oo, 1] can be associated fl], In particular, with 
n = 2cos(z^7r), by choosing u rational and fine-tuning the potential of the model one 
can reach any rational conformal matter field [^, |^. Until recently a solution of the 
model away from its critical points was known only for n = ±2 P, Q and for various 
rational values of z/ and potentials of low degree @, §]• In reference |^ we wrote 
down an exact solution of the model valid for any potential and any value of n. The 
solution was parametrized in terms of an auxiliary function. In the present paper 
we will determine this function explicitly as a combination of 6'-functions, thereby 
completing our solution of the model. Eventually the determination of the auxiliary 
function might in addition lead to a better understanding of the underlying continuum 
physics. In the one matrix model {n = 0) case the auxiliary function is the well known 
"Wigner square root" which was essential in revealing the Virasoro structure of the 
model and establishing its connection with the r-function of the kdV hierarchy ^ . In 
the general case information about the integrable structure underlying the 0{n) model 
should likewise be encoded in the auxihary function. 

Here we will in particular be concerned with using our exact solution to investigate 
so far unexplored regions of the parameter space of the model. Although our solution 
allows us to investigate any version of the model we will here restrict ourselves to the 
simplest "gaussian" one. In addition we shall consider only genus zero. The results 
can easily be extended to higher genera by means of the iterative procedure described 
in reference . We shall derive exact expressions for a number of quantities, including 
the eigenvalue distribution, away from the traditionally studied critical points. In 
particular we will use our solution to study the case \n\ > 2. As we shall see the model 
is well defined in a certain region of the coupling constant space both for n < —2 and 
for n > 2. For n < —2 we find no new critical points while for > 2 we find new 
critical points at which the string susceptibihty exponent •jstr takes the value +|, the 
value characteristic of branched polymers. 

We start by, in section 2, presenting the model and those of our results which will 
be of importance for what follows. We then proceed to, in section 3, determining ex- 
plicitly the above mentioned auxiliary function. In section 4 we outline the strategy for 
analysing completely any given version of the model and in section 5 we derive a closed 
expression for the second derivative of the free energy, an expression which proves very 
convenient when it comes to the investigation of the critical behaviour of the model. 
Hereafter we specialize to the simplest version of the model and section 6 contains a 



2 



detailed analysis of this version for all values of n g] — cxd, +cx3[. Finally, section 7 
contains our conclusion and a discussion of possible future directions of investigation. 



2 The Model 

The 0(?2)-model coupled to 2D gravity is described by the following partition func- 
tion m 



Z = e^'^ = [ dM ft dA, exp I -N Tr 

and Ai, i = 1, . . . , n are her 
potential to be completely general, i.e.Q 



n 

2 



(2.1) 



where M and A^, i = 1, . . . , n are hermitian N x N matrices. We shall take the 



V{M) = -M' (2.2) 
j=i J 

and use the following parametrization of n 

n = 2 cos(z/7r). (2.3) 

Let us summarize those of our previous results which are needed for the following. In 
the solution of the model the central object is the genus zero 1-loop function W{p) (or 
equivalently the distribution function p(A) for the eigenvalues {A} of the matrix M at 
the stationary point of the integral ( |2.1| )) 



W{p) = lim (^tr -^-rj) = I -^rfA. (2.4) 



A^-oo W p - M J p - X 

Once W{p) is known any other correlator of the M-field as well as the free energy as- 
sociated with surfaces of arbitrary topology can be found by an (in principle) straight- 
forward iterative procedure. In the following we consider the situation where the 
eigenvalues are restricted to only one interval [a, b] with a > and assume that the 
corresponding distribution function is normalized to 1. Under these circumstances the 
1-loop function W{p) is analytic in the complex plane except for one cut [a, b] and the 
eigenvalue distribution p(A) can be found from 

p{X) = ^{W{X-iO) -W{X + iO)} . (2.5) 
The 1-loop function fulfills the saddle point equation []1| 

W{p + iO) + W{p - iO) + nW{-p) = V'{p) (2.6) 

^We note that the traditional form of the partition function where there is no term hnear in M in 
the potential and a mass term for the A-field can be obtained by a linear shift of M. 
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which together with the requirement 

W(p)^-, p^oo (2.7) 
P 

determine uniquely W{p). The key idea which allowed us to solve ( |2.6| ) was a trans- 
formation of all functions involved in the problem; namely corresponding to a function 
h{p) we introduced (following reference 0) h±{p) by 

K{P) = 1^—^^ , h.{p) = h^{-p). (2.8) 

Then ( |2.6|) could be written in a more manageable form, W±{p) could be found and 
W{p) reconstructed using the inverse of (p78|) 

h{p) = -I (e*'^"/2/^+(p) - e-'^'^l'^h^ip)) . (2.9) 

The p-dependence of any multi-loop correlator of the M-field on surfaces of any genus 
could be described using a set of basis functions G'^^^ip) (and gI^\p)) defined by 

1. G'^^\p) and Gl^\p) satisfy the homogeneous saddle point equation 

G^''\p + iO) + G^''\p - iO) + nG^^\-p) = 0, p G [a, h] (2-10) 

2. G^^\p) and G'^^\p) behave near the end points of the cut [a, h] as 

Gi^)(p) ~ (p - ar^~"\p - b)-^'\ Gt\p) r^ip- b)-'^~y\p - a)-'/\ 

3. G^^\p) and g'^^\p) are analytical outside the cut (especially near —a and —h). 

4. G^^\p) and gI''\p) have the following asymptotic behaviour 

G'^:\p),Gi'\p)r.J^^, P-.00. 

To determine G'^^^p) it suffices to know G'f\p) = Gl^\p) = G'*^°^(p); namely using 
the important observation that any solution of the saddle point equation ( p.lO| ) can 
be parametrized in terms of any two other independent solutions, a recursive strategy 
for determining G'^^\p) from the knowledge of G'^'^\p) was formulated. Needless to 
say that G'^^\p) appears from G'^^\p) by the interchangement a ^ b. These basis 
functions contain no explicit reference to the matrix model coupling constants, {gi}- 
They depend on these only implicitly via the endpoints of the cut, a and b. The nature 
of this dependence is different for different values of n but this difference can be hidden 
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by working with a parameter, e, which is the only point (apart from oo) at which 
G^^\p) vanishes, i.e. 

(e) = GL°H-e) =0. (2.11) 

The exphcit dependence of the observables on the matrix model coupling constants 
could conveniently be described using a set of moment variables generalizing those 
introduced in reference [|1^ 



Mk= { p:V'{u)Gi'\u;), Jk = Mk{a ^ h) (2.12) 

where for n = 2 cos(z/7r) the G-f unctions are the G-functions corresponding to n = 
2cos((l — z/)7r) and where Ci is a curve which encircles the cut [a,b]. (We use the 
convention that all contours are oriented counter-clockwise.) As in the 1-matrix model 
case the moment variables have the advantage that for all observables of the model 
except the genus zero contribution to the one-loop correlator and the free energy the 
dependence on the infinite series of coupling constants {gi} arranges into a dependence 
on only a finite number of moment variables. Furthermore the moment variables reflect 
more directly than the coupling constants the possible critical behaviour of the model. 

In our previous paper we wrote down a contour integral representation for the 
genus zero contribution to the one-loop correlator, the integrand being composed of the 
functions G^'^\uj), G^^\uj) and the derivative of the potential V'{uj). For the following 
analysis it is convenient to rewrite our solution in the moment formulation. Doing so 
will also highlight the statement made above concerning the advantage of the moment 
variables. First let us remind the reader that the 1-loop correlator can naturally be 
decomposed into a regular part Wr{p) and a singular part Ws{p) as 

W{p) = Wr{p)-Ws{p) (2.13) 

with Wr{p) being given by 



4 — ra^ 



Now, let us introduce a scalar product by 

{f,g) = {n'-A)i ^{Uico)g^iu)-Uu)g4u;)} (2.15) 
Jc2 47r 

where G2 is a contour which encircles the two intervals [a, b] and [—b, —a] . Then it 
holds that (cf. to reference [0) 

Mk = {Ws,G^^^). (2.16) 
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Furthermore, let us define functions Wk{p) wliicli are solutions of the saddle point 
equation ( |2.10| ), analytic outside the cut [a,b] and fulfill 

{Wk,&I^) = 6,,i. (2.17) 

It is easy to see that Wk{p) must have the following asymptotic behaviour 

^kip)-^-. ^ j7z,p\ P^oo (2.18) 

4cos^ (?7fc7r/2) 

where 

V2k = i^, V2k+i = 1 - z/ (2.19) 
and that in the vicinity of a and b 

Wk{p)r^{p-a)'-'/\p-by/'. (2.20) 
This determines Wk{p) uniquely and one can verify that 

W,+ ip) = - ] A a'-p')'-'VpGt'\p) (2.21) 

^ S1H( Z-^TT ) 

where 



^= ^(p2_a2)(p2_62), (2.22) 

Here and in the following it is understood that we choose the sign of the square root 
so that y/p — s> oo as p — s> oo. It now follows that the singular part of the one- loop 



correlator can be written as 



Ws{p) = Y.MkWk{p) (2.23) 

k=l 



Mo = and M__i = I —u V'{uj) G^^\uj) = 2 - n. (2.24) 



and that the two boundary conditions which determine the endpoints of the cut a and 
b read 

ICi 

In particular from ( |2.21|) and ( |2.23| ) we can read off the eigenvalue distribution for any 
given value of u and any potential V{M). 

It is important to note that the decisive step in solving the model is the determina- 
tion of the auxiliary function G^^\p). Once this function is known all other quantities 
can be found. In our previous paper we derived a first order differential equation for 
G^^\p) and formally wrote down its solution. In the next section we will show, using 
a totally different strategy, how G^'^^ (p) can be completely explicited as a combination 
of ^-functions. This will make complete our solution of the model and explain sev- 
eral of our previous observations. In particular, deriving the explicit expression for 
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G^'^\p) gives as a byproduct a determination of the parameter e in terms of z/, a and 
b. Furthermore the observation that when u is rational G^^\p) reduces to an algebraic 
function is an immediate consequence of the general formulas. 

Before proceeding with the determination of G^^\p), let us spend a few lines showing 
how the possible types of critical behaviour are very clearly exposed in the moment 
formulation. For u G [0, 1] critical points are reached when the endpoint a of the 
eigenvalue distribution touches the origin. In this limit one has (cf. to reference [0) 

^p-'^-'?^, p^O (2.25) 

which means that 

W^fc(p) (2.26) 
Hence we have a series of M'th multi-critical points characterized by a critical exponent 

Ws (p) ~ p''^' , /XAf = M - rim+1 , (2.27) 
with the corresponding subspace of the coupling constant space being given by 

Mo = Ml = . . . = M„_i = 0, M„ ^ 0. (2.28) 

This reproduces the well-known results of references |jl|, |^. 

3 Determination of the function G^^\p) 
3.1 Reformulation of the problem 

To determine the function G^'^\p) it is convenient to perform a change of variables 
from p to u given by 

p = asia{u,k), k = -. (3.1) 

This maps the complex p-plane into the rectangle [—K^ K] x [—iK', iK'] in the complex 
M-plane. An important feature of the mapping ( p.l| ) is that it "opens up" the cut [a, b] 
of G^^\p), i.e. the part of the complex p-plane which lies above the cut [a, b] is mapped 
into the first quadrant of the M-plane while the part of the complex p-plane which 
lies beneath the cut [a, b] is mapped into the fourth quadrant of the w-plane. The 
analyticity properties of G^^^ (p) imply that the function G^^^ (u) must fulfill a number 
of relations on the boundary of its domain. First of all, the saddle point equation ( |2.10| ) 
implies that for u G [—iK', iK'] we have 

G^^\K + u) + G^'^^K -u)+ nG^''\-K + u) = 0. (3.2) 
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Furthermore the fact that G^'^\p) as no cut along the interval [—6, —a] gives that for 
u G [-iK',iK'] 

G^'^^-K + u) + G^^\-K -u) = 0. (3.3) 

Finally from the fact that G^^\p) is analytic along the line segments [—oo,—b] and 
[b, oo] we deduce that for u G [—K, K] 

G'-^^tK' + u) = G^''\~iK' + u). (3.4) 

By means of the three relations ( |3.2| ), ( |3.3| ) and ( p.4| ) we now analytically continue the 
function G'('^^(u) from its original domain [—K,K] x [—iK',iK'] to the whole complex 
plane. First we extend the definition of the function G('^^(u) to the vertical band 
[—K, K] X [—zoo, ioo] by means of the relation (|3.4|). Next we use the relation ( |3.3| ) to 
extend the definition of the function to the band [SK, K] x [—ioo, ioo] and then using 
the equation ( ^.2[ ) we can define G^^\u) in the remaining part of the complex plane. 
This procedure leaves us with a function which is defined in the whole complex plane 



and which obeys the equations ( p. 21) , (|3.3| ) and (|3.4|) for all u. Now we would like to 
determine this function. First, let us note that combining (p.2|) and (p.3|) we get 



{u + 2K) + {u - 2K) + nG^^^ {u) = 0. (3.5) 



Furthermore the parity condition (|3.3|) and the periodicity condition (|3.4| ) can be ex- 
pressed as 

G(°)(-n-2ir) = GW(«), (3-6) 

G^''\u + 2tK') = G(°)(m). (3.7) 

Let us next introduce a translation operator X by 

X = e2^^. (3.8) 

Then we can write ( p.5|) as 

(X^ + nX + l) G(°)(w) = - {X - x_) G(°)(w) = (3.9) 

with 

= eT^^i-'^)-. (3.10) 
Now, let us decompose G^'^\u) in the following way. 

G^^\u) = ^Gf{u) + y/^G^°\u) (3.11) 
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with 



Gf(u) = — (X -X-)G^^\u), (3.12) 



(0) 



u) = — i (X -x+)G^^\u). (3.13) 



G 

Then it obviously holds that 

XGf\u) = x+Gf\u), XG^^\u) = x^G'^°\u) (3.14) 

and it is easy to show that the decomposition ( p.ll|) of G^^\u) is unique, i.e. if G^^^ = 
iJ+ + with XH± = x±H± then H± = yyx±G±\ Exploiting this uniqueness result 
it follows from the parity condition ( p.3|) that 

gL°Hm) =Gf (-«)• (3.15) 

It is now clear that the decomposition ( ^.ll] ) corresponds exactly to the decomposition 
introduced in equation ( p.9| ) and our original problem of determining G^^^ (p) has hence 
been transformed into the problem of determining a function ■* (u) which fulfills 

Gf\u + 2K) = e-'^^-'^^Gf{u), (3.16) 
G^+\u + 2iK') = G^+\u) (3.17) 



(where the second equation originates from the periodicity condition ( |3.4| )) and which 
is compatible with the requirements 2, 3 and 4 on page ^. In the following subsection 
we will show how to solve this problem. 

3.2 The explicit expression for G^^\p) 

Let us start by noting that the condition 2 and 3 on page | imply that G^^\p) has 
simple poles at u = K and u = K + iK' but not other singularities. Furthermore the 
condition 4 implies that G^^\u) has a simple zero at u = iK'. Now we introduce a 
function H^{u) by 

Q f U — iK' \ Q ( U — £ 



H+{u) = ^-^T^^^-^^^^e-^i-'^)^. (3.18) 



-{K+iK') 
2K 



This function has the same poles as (?+''(«). It also has the zero at m = iK' in common 
with G^^{u) and it fulfills the equations (|3.16|) and (|3.17| ) provided the parameter e 
takes the following value 

e = i{l-v)K'. (3.19) 
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Now, Gf{u)/H+ (u) is a doubly periodic function (with periods 2K and 2iK') which 
has no singularities. (The function can have no singularities except for poles and of 
poles it can have only one originating from the zero u = e of H^{u). Hence it has none 
at all.) This allows us to conclude that 

G^+\u) = const ■ H+{u). (3.20) 

In particular we see that in accordance with our previous analysis the function (u) 
must have exactly one zero in addition to the one at m = iK' {p = oo), namely u = e. 
Translating back to p- variables gives us the value of e (cf. to relation (|2.11|) ) 

e = asn{i{l-u)K'). (3.21) 



Now it only remains to determine the constant in relation (|3.20|) . Its value follows from 
the normalization condition 4 on page ^ and reads 

const = y — ^ — . (3.22) 

a ^4^4(2^) 2cos(i^7r/2) 

Now we have the explicit expression for G^_^\u) and we can write down as explicit 
expressions for (in principle) any other quantity. Let us note in this connection that 
with our sign convention for the square root the quantity translates to m- variables 
as 

= — a6 cn-u dnti (3.23) 

3.3 Rational case 

Let us remind the reader that when u is rational the scaling behaviour at the critical 
points of the model is that characteristic of 2D gravity interacting with rational con- 
formal matter fields. More precisely, when u = l/q with I and q integer and < / < g 
the matter fields which appear are of the type (g, (2m + l)q±l). In reference it was 
shown that when u is rational a simplification of the function G^^\p) occurs; namely 
foT u = l/q one has 

{VpG^ip))' = Mp)Vp + B{p) (3.24) 
where A{p) and B{p) are polynomials in p having a definite parity 

A{-p) = i-iy^'Aip), B{-p) = {-iyB{p) (3.25) 

and having degree (g — 3) and q when {q + I) is even and {q — 2) and {q — 1) when 
(g + Z) is odd. The existence of a relation like (p.24|) is a direct consequence of standard 
properties of elliptic functions. Consider the identity ( p.24| ) rewritten in the w-variable 

(cnu dnu G'"^\u)y = A{snu)cnudnu + B{snu). (3.26) 
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From ( |3.16| ) and ( |3.17| ) it follows that the function on the left hand side of ( ^.261 ) is an 



elliptic function with periods 2K and 2iK' when I is even and periods AK and 2iK' 
when I is odd. The same is true for the function on the right hand side. Furthermore 
the two functions have the same poles, namely only one of order q ai u = iK' Now, 
counting the number of adjustable constants in the polynomials A and B it is easily seen 
that one can always arrange that the two functions also have the same zeros namely one 
of order q aX p = e. Exact equality between the two functions can hereafter be ensured 
by choosing the overall normalization of A and B so that the residue at m = iK' equals 
1. 

4 Explicit solution of a given model 

Our explicit expression for the auxiliary function G'^'^\p) gives us the possibility of 
exploring in detail the coupling constant space of our model. In particular exact results 
for a large number of quantities can be obtained. Let us briefly describe how one would 
extract exact results for a given model. Observables generally depend explicitly on the 
matrix model coupling constants {gi} via the moments and implicitly via the endpoints 
of the cut a and h which are in turn determined by the boundary conditions ( p.24| ). 
For a given matrix model potential the moments (including those entering ( p.24| )) 
can be expressed in terms of the a and h by simply carrying out the contour 

integrals appearing in the definition (|2.12 ). This can either be done using the 9- 



function representation of the G-functions in which case the contour integration can 
be reduced to to an integration around the pole u = iK' or using the p-representation 
in which case the contour can be deformed into one which encircles infinity. We will 
choose the latter line of action. To proceed along this line one needs to know the large 
p-expansion of G^^\p) 

(^(o)fp\-^fl ^tan(z/7r/2)ai 0^2 i tan (z/7r/2) gg \ 
p \ p p^ p^ ) 

As explained earlier once G'^'^\p) is known the remaining G-functions can be found by 
a straightforward iterative procedure. The a-coefficients are most easily found from 
the differential equation for G^^{p) derived in reference 



^ [VpGT{p)) = I «i - — I G'^'^\p)+pgMGV{p) (4.2) 



where 

g+{p) = - 

p^ — e 



gM = ^^t^f • (4-3) 
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Here the constant ai is an integration constant which is not determined by the differ- 
ential equation itself but can be found by other means It is given by 



e^/e de^ a? {I? — 
and fulfills the following relation which we shall make use of later 



In the following we will use the notation that the expansion coefficients of the function 
&^\p) are denoted as {cui}. Obviously the {ai} appear from the {ai} by the replace- 
ment u 1 — u. The relationship between the two sets of parameters can also be 
expressed in another often very useful way, namely via the identity (cf. to reference ^) 

G^^\p) = - cot (^) g+{p)Gf{p). (4.6) 

For instance this identity implies 

a, = (^a^ - ^ j . (4.7) 

Let us in this connection also note that we have the relations (with the obvious nota- 
tion) 

ab ^/e \/I 



(4.8) 

e e e 

as well as 

Pa = l- pa= pb- (4.9) 

Here the first two relations are simple consequences of the identity (|3.21| ). (We note that 
the sign of -y/e is determined by equation ( |3.23| )). The first equality sign in the second 
relation follows from the relation e = —ab/e while the second follows from the fact that 
the parameter ai should be symmetric in a and b. Inserting the expression (|3.21|) for 
e in ( [4. 41 ) we get for ^ 

^ = E{ivK') + 'w{E' - K') = Z{wK') + iu^. (4.10) 
b 2K 

Furthermore we find for pa 

For our later consideration we shall also need to know a2- It reads 



a, 2 



-(e^-al-a^-b''). (4.12) 
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In principle the boundary equations can be explicitly solved and the moments deter- 
mined for any matrix model potential but of course the boundary equations become 
more and more involved as the degree of the potential increases. Let us point out that 
our expressions are valid not only for any potential but also, at least formally, for any 
value of V. Hence our formulas can be used to explore the hitherto unexplored regions 
of the parameter space, n < —2 and n > 2. While a priori nothing prevents us from 
choosing these unconventional ranges for n there is of course no guarantee that results 
obtained for \n\ > 2 are physically meaningful. For example it would not be acceptable 
(at least immediately) if the endpoints of the cut turned out to be complex or if the 
eigenvalue distribution were not real and positive. 



5 The string susceptibility 

In this section we will derive a closed expression for the string susceptibility which we 
will make use of later when investigating the critical behaviour of the model. First let us 
introduce an overall coupling constant in front of our potential, i.e. let us replace V{p) 
by V{p)/T where T can be thought of as the cosmological constant or the temperature. 
The string susceptibility, U(T), is then given as U(T) = ^(T^Fq) where Fq is the 
genus zero contribution to the free energy defined in equation (|2.1| ). Knowing the string 
susceptibility we can, in case the model has a critical point at T = T^, extract the value 
of the critical index jstr 

U{T) = ^^{T'Fo)r^{n-Tr'- (5.1) 
To proceed with the derivation we first note that 



where the subscript refers to the genus zero contribution. Next we use the identity 
(proven in reference 0) 

-^{TWip)) = G^'\p) (5.3) 

to conclude that 

-^{t'Fo)=I^^^G^'\pMp) (5.4) 

and consequently that 

f^,dFo \ ^ ^dim ^ r dp ^^^^ r dG('\p) da' ^ dG^d^) 



dT' \ dT dT Jc, 2m '"^^ I da' dT db' dT 
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Let us rewrite ^'^g^i^^ in the following way 

^^ = |;M'°'(p)+<=2pG<««(rt). (5.6) 

The existence of such a relation follows from the fact that any solution of the homo- 
geneous saddle point equation (|2.1CI| ) can be parametrized in terms of any two other 



independent solutions (cf. to reference 0). The two constants ci and C2 can be deter- 
mined using the fact that the right hand side of equation ( |5.6| ) should have the same 
asymptotic behaviour and analyticity structure as the left hand side. The asymptotic 
behaviour is determined by (|4.1| ) and (|4.4|) and as regards the analyticity structure the 
requirement to be imposed is a behaviour of the type (p — 6)^^/^ as p b. We note 
that the determination of Ci and C2 is facilitated be the use of the relation (f4.6|) . The 
two constants read 

1 1 e^-a^ 

Now we can actually evaluate the integral (|5.5| ); namely inserting the expression (|5.6|) 
for ^'^g^2^^ and its equivalent for '^^gj^J^^ in (|5.5| ), integrating by parts and making use 
of the boundary equations we find 

dU{T)_(^ n\ 1 [e^-a^da^ - 6^ 1 

1-77 15 -5\ — 15 — -1^ U — -J^)- (5.9) 



dT \ 2) h'^-a'^ \ a? dT 6^ ^ 

This quantity is amazingly universal. First, it contains no explicit reference to the 
matrix model coupling constants {gi}- Secondly, the given expression is valid for all 
values of n. This universality is another manifestation of the universality of the two- 
loop correlator of the model observed earlier |]^. For = the universality with 
respect to the coupling constants was discovered in reference |jTl|. In that case the 
above equation can be integrated exactly which gives 

t/(r) = 21og(6-a), n = 0. (5.10) 

We note that for n = 2 special care must be taken. The parameter e diverges and 
the prefactor (^1 — goes to zero. The case n = 2 is not by any means a singular 
case, however. Only our parametrization is not well suited for this value of n. Actually 
as shown in references P, |, |^ for n = ±2 a simpler parametrization can be chosen. 
The limit n — > 2 of all expressions in the present parametrization are well defined and 
reproduce the results obtained in the simpler parametrization. 
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6 The gaussian potential 



We will now solve in detail the 0{n) model in the case of a gaussian potential V{M). 
Due to the interaction term in the action this is by no means a trivial case. On the 
contrary one expects that the major features of the general model are reflected already 
in this simplest version. We choose our potential to be of the form 

Vip) = ^{p-po)\ Po>0 (6.1) 

and we remind the reader that in addition to the potential V{p) the eigenvalues of 
the matrix M feel an attractive {n > 0) or repulsive {n < 0) force from their mirror 
images with respect to zero and that when one of the eigenvalues touches the origin the 
partition function ( p.l| ) ceases to exist (cf. to reference [Q). With a potential like ( |6.1| ) 
one has a well around the point p = po and one expects that by choosing T sufficiently 
small one can obtain a stable situation where the eigenvalues are confined to the well 
and located at a finite distance from the origin. Stated differently one expects that a 
stable 1-cut solution of the matrix model exists in some region of the coupling constant 
space. In the following we shall explore in detail the coupling constant space of the 
model and determine when such a stable solution exists. In particular we shall look for 
singular points, i.e. points where the solution ceases to exist or changes its nature. In 
order for a solution to make sense the endpoints of the cut a and b which are determined 
by the boundary equations ( p.24| ) should come out real and positive. Let us take a look 
at these boundary equations. Deforming the contours in ( p. 241) to infinity and inserting 
the large p expansion for G^^^ (p) and G^^-* (p) we find that the boundary equations with 
the potential ( |6.1| ) read 

ai + ipo tan (^^^ = 0, (6.2) 



p'q tan' i^—j +e' = 2(2 - n)T (6.3) 

where we on our way have made use of the relations ( [4.7D and ( [4.12|) . For given values 
of the parameters of the potential these equations determine a and b. However, due 
to the complexity of the equations (cf. to ( |3.21| ) and ( [4.10| )) trying to solve directly 
for a{po,T) and b{po,T) is not a practicable way of proceeding. We shall hence take 
another line of action. Inserting the expressions ( [4.7[ ) and ( ^.21| ) for ai and e in ( |6.2| ) 
and ( |6.3| ) and performing a few additional manipulations we arrive at the following two 
relations 

ipo tan ( ^ ) 

b = — (6 4) 

E{iuK') + iu{E' - K') ' ^ ' 
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Let us fix the parameter po- This corresponds to fixing the position of the potential well 
and does not influence qualitatively the features of the model. Then we have that for a 
given value of the parameter = f the quantities b, T and a are uniquely determined. 
Our strategy for studying the above equations will hence be the following. We let k 
vary between and 1 and determine for each of its values the corresponding values of 
a, b and T. In doing so we find a(T) and 6(T) and miss no real solutions. It is easy to 
show that the relations (|6.4|) and ( |6.5|) give real values for a, b and T for any k G [0, 1] 
and any n g] — cx3,oo[. These quantities are in addition positive except for certain 
ranges of k values when g] — oo, — 2[. We shall explain the situation in more detail 
later when we consider separately various ranges of n. Let us just note that bearing in 
mind the parametrization of G'-"^ (p) is section |^ the use of the quantity = | as the 
fundamental parameter is actually quite natural. 

Let us now turn to determining the eigenvalue distribution. From ( ^.2tj| ) and ( |2.12| ) 
it follows that we have 

W,{p) = ^W,{p) (6.6) 



and hence according to (p.5|) and ( |2.13|) 



p(A) = ^W,iX) 



X+iO 



X-iO 



(6.7) 



Using (U), ( ^^211) and the rotated version of the saddle point equation ( p.lO|) (cf. to 



reference 0) we can rewrite ( |6?7| ) in the following way 

Now bearing in mind the properties of the mapping (|3.1|) this identity can also be 
expressed as 



A+iO 

(6.8) 

A-iO 



• _ K—iv 

p{v{\)) = —— e*^"/2 ^^5) (m) 

2nT K+iv 



(6.9) 



where v is related to A by 

A = asn(K + if), ve[Q,K'] (6.10) 

Finally if we insert the explicit expression for (5+^ {u) found in section ^ we arrive at 
the following expression for the eigenvalue distribution 



2vrT^ 2sin(f)^4(W) 



^2 e - 6»2 — — e^^K \ (6.11) 



^3(^)1 n J '\ 2K 
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It is worthwhile noting that the expression ( p.ll| ) is vahd for any value of v. In particu- 



lar ( |6.11| ) together with the two boundary equations ( |6.2| ) and ( |6.3| ) give a generalization 
of the Wigner semi-circle law which is reproduced when n is set equal to zero (z/ = 1/2). 
As mentioned earlier we can of course not be sure that the solution makes sense for 
any range of the parameters v and T. We shall address this aspect in detail later. 

Let us finish this section by writing down the expressions for da^/dT and db'^/dT 
for the gaussian potential ( |6Jl) . These can be found by differentiating (|6^) and (|6.3|) 
with respect to T and making use of the relation (|4.5| ) and ( [4.9| ). They read 

(6.12) 



da^ 2(2 -n)a^-e^ 1 db^ 2(2 - n) b^ - 1 



dT T a^-b^l-pa dT T b^ - 

which implies that ^^^^ takes the form 
dU{T) (2 - nf 



dT T (62 




(6.13) 



From this expression we can readily read off where in the coupling constant space we 
can expect to encounter critical points; namely we see that the free energy can become 
singular when pa is equal to zero or one, when it diverges, when ^ oo and due to 
the relation ( p.21|) when a — 0. 



6.1 nG]-2,2[ 

For n g] — 2, 2[ we have < u < 1 and we can rewrite the equations (|6.4| ) and ( |6.5|) as 

b = - (6 14) 

uE' - E{uK', k') + dn(z/is:', k') sc{uK', k')' ^ ' ' 

T = ^^^{p^an(^)-A:Vsc^(.ir',fc')}. (6.15) 



It is straightforward to show that the expressions ( p.l4| ) and ( |6.15| ) give real non nega- 



tive values for b, T and a {= k-b) for all k G [0, 1]. In figure 1 we have plotted a{T) and 
b{T) for various values of u. The curves are obtained by letting k vary between zero 
and one. We have set po equal to one. Note that the point a = b = po corresponds to 
k = 1 and the point a = to A; = 0. In accordance with the nature of the interaction 
between an eigenvalue and its mirror image, we see that the values of a and b increase 
with u for fixed T. (This can also easily be inferred from the boundary equations.) Let 
us stress that for any given value of k we have a closed expression for the eigenvalue 
distribution, namely (|6.11| ). 

Let us show how our formalism allows us to recover immediately the results of 
references 0, ^, ^ ^ regarding the critical behaviour of the model for the here 
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T 



Figure 1: The variation of a (lower curves) and h (upper curves) as a function of T for 
V = 1/10 (full line), v = 1/2 (dashed line) and z/ = 2/3 (dotted line). 



considered range of n-values. As mentioned earlier the free energy can become singular 
if pa becomes equal to zero or one, if it diverges, if ^ oo or if a ^ 0. From the 
parametrization ( p.21|) it is obvious that remains finite for v €]0,1[ and from the 
expression ( ^4.11| ) for pa it follows that 

< Pa < 1, for 0<v <l. (6.16) 

(Actually Pa = 1 for A; = 1 for all values of v g]0, 1[ but this case is not interesting from 
the point of view of critical behaviour since it corresponds to a = h = p^^T = 0.) Hence 
the only possible type of singular behaviour is associated with a = 0. Using ( |6.14| ) 
and ( |6.15| ) we find the corresponding critical values of h and T 

tan f ^) 1 
z/ 2(2 + n) 



From equation (|3.21|) it follows that in the limit a — > 0, e behaves as 



a 



e^2ih[^-j (6.18) 
which according to (|0| ) implies that 

U{T) ~ a^^ (6.19) 
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Now, since Pa ^ I-' (cf. to equation ( |6.18| )) we get using ( |6.12| ) 



(6.20) 



Note that this scahng law is very clearly exposed in figure 1. The relation (|6.20|) 
together with ( |6.19D allow us to extract the value of 7^4^ associated with the critical 
point a = 0. We find 

Istr = (6-21) 



which indeed coincides with the result of references |jl|, ||, |^ . In addition, using the 
expression ( |6.11| ) one can also easily recover the expression for the eigenvalue density 
at the critical point found in reference 0. 



6.2 



n 



±2 



For n = ±2 the equations ( |6.14| ) and ( |6.15| ) contain divergent factors and special care 
must be taken. As mentioned earlier for these values of n the present parametrization 
is not the optimal one and a simpler one can be chosen. However, the limits n — > ±2 
of ( |6.14| ) and ( |6.15| ) are well defined and lead to the same equations as are obtained in 
the simpler parametrization j^, |l| 

6.2.1 n = -2 

For n = —2 we find from ( |6.14|) and (|6.15| ) taking the appropriate limit 



^X', T = - de + 1)K' 



2E' 



(6.22) 



Furthermore we see from (|3.21|) that the parameter e is equal to zero and hence we can 
integrate exactly the expression in equation (|5.9|) which gives 



(6.23) 



U{T) = 2log{b^ -a^). 



This is in accordance with the observation that the solution of the n = —2 case can 
be read off from the solution of the n = case (cf. to references @, 0)- From ( |6.23| ) 
and ( |6.12| ) it follows that singular behaviour can only occur if pa becomes equal to zero 
or one. For pa we find using ([4.11|) 



pa 



k 



E' 



(6.24) 



This quantity is always less than one but becomes equal to zero as /c ^ 0. However, 
as is seen from ( |6.22| ) the point k = corresponds to 6 = T = 00 and is not interesting 
from the point of view of continuum behaviour. We note that be = = 00 is also what 
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one would expect from the formulas ( |6.17D . We also note that it is due to our choice 
of a gaussian potential that the model has no critical point. It is well known that for 
potentials of degree larger than two there exist critical points which are characterized 
by having logarithmic scaling violations 

6.2.2 n = 1 

For n = 2 the expressions (|6.14| ) and ( |6.15|) turn into 



2E^ 



(6.25) 



Furthermore, taking the limit n ^ 2 of equation ( |6.13D leaves one with an expression 

IU{T 
dT 



for which can be integrated exactly and leads to the following result 



UiT) = 2(^] log 



(6.26) 



\2K'J '"^\b'^~a\ 

Obviously a = is a singular point. The corresponding critical values of b and T read 



(6.27) 



which is again in compliance with the general results (|6.17|) . For the quantity pa we 
find using ( |4.11D 

1 r F'l 

(6.28) 



E' 



It is easy to see that < pa < 1/2 and that = iff /c = 0. Taking the limit a — 
we find 



U{T) 



where according to ( |6.25| ) 



log a 

T,-T ^ a\\ogaf 



(6.29) 
(6.30) 



6.3 n > 2 

For n > 2 we sei V = iv with u real. Hence we have 

n = 2 cosh(z/7r). 

With this parametrization we can write our two boundary equations 

Potanh(f ) 
^ " E{pK') + u{E' - K') ' 

^ - ^{pSta„h(f)-.W(.A-')}, 



(6.31) 
and ( |6.5D as 

(6.32) 
(6.33) 
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Figure 2: The parametric curves {T{k), a{k)) and {T{k), b{k)) for k g]0, 1] and u = 1.2. 



Also in this case one can easily convince oneself that the expressions ( |6.32| ) and ( |6.33| ) 
give rise to real non negative values of b and T for all k G [0, 1]. Furthermore for each 
value of 9 there is a maximal temperature, namely 

pI (6.34) 



2(2 + n) 



We note the similarity with the expression for the critical value of T, ( |6.17| ) for n G 
] — 2, 2[. However, in the present case the maximum temperature in attained not only 
once but an infinite number of times corresponding to the infinite number of k values 
which solve the equations 

uK' = 2mK, m = l,2, .... (6.35) 

In other words for k larger than some critical value kc given by 

K' 2 

-(fc,) = (6.36) 
A u 

an oscillation of the temperature begins. In figure 2 we have shown a{T) and b{T) for 
u = 1.2. The curves have been produced by letting k vary between zero and one and 
as before po has been set equal to one. In the parametric plot we see an oscillation of 
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Figure 3: The parametric curves {T{k),b{k)) and {T{k),a{k)) for k G [0.02; 0.38] and 
z/ = 1.2. 



the endpoints a and 6 as a function of T. (Due to finite precision only a finite number 
of oscillations appear.) The right turning points for the oscillation correspond to the 
/c- values which fulfill ( |6.35D . At these points both ^ and ^ diverge. This divergence 



is due to the divergence of e. (We note that pa also diverges at these points.) The 
left turning points correspond to the situation pa = 1. The associated k values are the 
solutions to the equation 

E{^K') + u{E' - K') = ^ ^^^_^)^ ^ . (6.37) 

At these turning points ^ diverges but ^ stays finite. The value of T is the same 
at all right turning points but the value of T at the left turning points increases as k 

dT 

— n TVioir 

dT 

vK' = {2m + l)K. (6.38) 



decreases. We note that there is no divergence of ^ when pa tends to zero. On the 
contrary these are the points where ^ = 0. They are given by 



The vanishing of ^ at these points is due to the fact that e becomes equal to b. There 
is no singular behaviour associated with these points. In figure 3 we have enlarged 
the region of figure 2 where the oscillations of a(T) and 6(T) occur and we see that 
the features of the curves are in accordance with our analysis. This picture of course 
implies that at a certain temperature, Tq, namely the temperature at the first of the 
left turning points, the boundary equations start to have more than one solution and 
as T approaches Tmax the number of possible solutions goes to infinity. One may 
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wonder if there is a criterion that would allow one to pick out one of these solutions. 
We will argue that there is; we conjecture that only the solutions corresponding to 
the upper branch of the curves for a(T) and b{T) give rise to a positive eigenvalue 
distribution. In other words the model is only well defined if < k < 1 where kc is 
given by ( |6.36| ). To support this statement, let us determine the eigenvalue distribution 
at various selected values of k. First, let us consider the k values given by ( |6.38| ). As 



mentioned above at these values of k one has ^ = 0. Obviously the first such k gives 
a set of points (a(T),6(T)) which lie on the upper branches of the curves for a(T) 
and b{T). In general one will encounter one of the points corresponding to (|6.38|) each 



time one moves from a left turning point to a right turning point in the direction of 
decreasing k (cf. to figure 3). Inserting the equality ( |6.38| ) in the expression (|6.11| ) for 
the eigenvalue distribution we get 

, ,,,, 1 , .,rr, 1 / TIV \ sn(v,k') 

p(v(X))= ai-l)"^ cos u— , \ ,{ 6.39 

^ ttT ^ ' sinh(f) \2Kjdniv,k') ^ ^ 

where A is related to v by 

X = asn{K + iv), ve[0,K']. (6.40) 



Bearing in mind the relation ( |6.38D it is obvious that the eigenvalue distribution given 
by (|6.39|) remains positive for all v in the given interval only if m = i.e. only for the 
value of k which corresponds to a point on the upper branch of the solution. Next, 
let us determine the eigenvalue distribution corresponding to the right turning points 
given by ( |6.35| ). Inserting (|6.35|) in ( p.ll| ) gives 



p (v(X)) = 4^ a (-1^+^ sin ( jy^) ^ / , , (6.41) 

^ ttT ^ ^ sinhff) \2Kjdniv,k') ^ ^ 



with V being related to A as in ( p.40| ). In this case we see that the eigenvalue dis- 
tribution is positive only if m = 1 i.e. only at the first of the right turning points. 
Furthermore ( |6.39| ) and ( |6.41| ) show that as k decreases the eigenvalue distribution 
oscillates more and more rapidly from positive to negative values, the number of os- 
cillations going to infinity as A; — 0. In particular it is obvious that the model makes 
no sense at the point k = which would be the naive analytical continuation of the 
critical point for n e] — 2,2[to n > 2. 

However, now that we have rendered probable that the model is actually well defined 
for kc < k < 1 we shall determine the critical index jstr associated with the obviously 
singular point k = kc- For that purpose we let k —* {kc)+ and denote the value of T 
corresponding to kc as Tc. (We note that Tc = T^ax (cf. to equation (|6.34| )). From the 
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boundary equation (|0|) we see that 



e2 



Then introducing e in stead of e in relation ( p.l3| ) we get 

Now, from the relation (|4.7| ) we find 

ib' -a')il-pa) = ^a + b'-~e'^c- (T, - T)-'/\ (6.44) 
e 

We note that the constant, c, entering this relation is positive as it should be. The 



quantity a is negative as can be seen from the relation (|6.4| ). The parameter e is 
likewise negative which follows from the parametrization ( p.21| ) and the fact that we 
study the limit k — > {kc)+ with kc given by (|6.36| ). Finally is a positive quantity 



which is a consequence of the recipe (|]2|) and the identity (^M )- (We note that 
the sign of a/I changes from one right turning point to the next.) Combining ( |6.44| ) 
and ( |6.43| ) we get 

^~(T.-T)"3/^ (6.45) 

which means that 

Istr = +^ (6.46) 

6.4 n< -2 

For n < —2 we set u = 1 + iu with u real. Hence we have 

n = -2 cosh(i>7r) (6.47) 
and we can write the two boundary conditions ( |6.4| ) and ( |6.5D as 

po cothff ) 

E{uK') + cn{uK') ds{uK') + u{E' - K'Y ^ ' ' 

Obviously these equations give real values for a, h and T for all values of /c G [0, 1]. In 
order to describe in more detail the behaviour of a, h and T as a function of k, let us 
take a look at some selected points. First, let us consider the k values characterized by 

vK' = (2m + 1)K, m = 0, 1, . . . (6.50) 
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These points are analogous to the points given by ( |6.38| ) for n > 2. Only the roles of 
a and b are interchanged. In the present case we have that = a^, 1 — = and we 
find that ^ vanishes while ^ stays finite. As in the case n > 2 there is no singular 
behaviour associated with these points. We note that a, b and T are always positive 
and finite when k fulfills (|6.50|) . Next, let us consider the points given by 

i>K' = 2mK, m = l,2,... (6.51) 

When one approaches such a point in the direction of decreasing (increasing) k one 
finds b 0_, T 0+, (6 0+, T — > 0+). Thus b changes sign at at these points. 
Obviously b must also change sign at a series of k values each of which lies in between 
two successive solutions to (|6.50|) and ( |6.51D . These k values are characterized by the 



denominator in the expression (|6.48| ) becoming zero, i.e. 

E{i)K') + cn{i)K') ds{i)K') + V{E' - K') = 0. (6.52) 

At these points Pa = and both a, b and T diverge. It is easy to convince oneself that 
T always stays positive and that there are no other points than those given by ( |6.51| ) 
and (|6.52|) where b can change sign. In summary, b starts out positive and equal to 
Pq at k = 1, goes to +oo as k approaches the first k which solves (|6.52| ). In the same 



interval T increases from to oo. At the first solution to ( |6.52| ) a cycle starts. When k 



varies between one solution to ( |6.52| ) and the next, b and T behave in the following way. 
Precisely at the solution to ( |6.52|) b jumps from +oo to — oo. Hereafter it increases 
as k decreases and tends to zero when k approaches a solution to ( |6.51| ). In the same 



/c-interval T decreases from oo to 0. Then b and T increase again and both go to +oo 
when k approaches the next solution to (|6.52|) . These cycles are the analogues of the 
oscillations of b{T) and a(T) seen for n > 2. In particular we are lead to the conclusion 
that as in the case n > 2, the model makes no sense in the point k = which would 
be the naive analytical continuation of the critical point for n G [—2, 2] to n < —2. 
However, unlike what was the situation for n > 2, in the present case there are no 
interesting new singular points. The points 1 — = are, as mentioned above, not 
singular and at the points p^ = both a, b and T diverge. The final possibility e oo, 
(pa — oo) corresponds to the points given by ( |6.51[ ) and here both a, b and T vanish. 



Even though we have now made clear that our model has no interesting new critical 
points for n < —2, let us spend a few lines discussing where in the coupling constant 
space the model has a meaning for this range of n's. From the analysis above it 
follows that for each value of T there exists an infinite number of solutions for a and 
b (even when we reject the obviously unphysical solutions with a and b negative). As 
in the case n > 2 we will argue that only the solutions corresponding to the first 
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branch of the curves a{T) and b{T) make sense, i.e. the model is only well defined for 
kmin < < 1 where kmin is the solution of the first of the equations ( |6.52| ). To support 
this statement, let us as before determine the eigenvalue distribution at a set of selected 
k values; namely the k values which solve the equations ( |b.5(J ). All of these give a and 



b but only the first one belongs to the interval [kmm, !]• Inserting the identity (|6.50|) in 
the expression (|6.11| ) for the eigenvalue density we find 



^ ttT^ ^ cosh(f)^ ^ \ 2KJ dniv,k') ^ > 

with V related to A by 

X = asn{K + w), ve[Q,K'] (6.54) 

In accordance with our statement the eigenvalue distribution is positive only at the 
first of the here considered k values. 

7 Conclusion 

Having determined explicitly the auxiliary function we have completed our solution 
of the 0{n) model on a random lattice. Our solution allows an exhaustive analysis 
of the model for any value of n and any potential. We have carried out this analysis 
for the simplest "gaussian" version of the model which corresponds to a collection of 
self-avoiding non-intersecting loops densely packed on a randomly triangulated surface. 
Our analysis showed that the model is well defined in a certain region of the coupling 
constant space both for n < —2 and n > 2. For n < —2 we found no new critical points 
while for n > 2 we found new critical points at which the string susceptibility exponent, 
7sfr takes the value +|. We expect to encounter the same situation if we include higher 
order terms in the potential, i.e. we expect that the model will still be well defined in 
a certain region of the coupling constant space and that there will be new, possibly 
multi-critical points for n > 2. It is tempting to speculate about a connection between 
n > 2 and 2D quantum gravity coupled to matter fields with c > 1. Unfortunately 
there does not exist any mapping of the 0{n) model for n > 2 onto a model which is 
known to have c > 1 on a regular lattice. Let us note anyway that recent numerical 
simulations show that ''jstr changes rather rapidly from to +| when one crosses the 
c = 1 barrier . 



A less speculative unclarified point concerns the relation of the 0{n) model on a 
random lattice with the theory of integrable hierarchies. As mentioned earlier, in the 
one matrix model (n = 0) case the function that we have denoted as our auxiliary 
function played an essential role in revealing the Virasoro structure of the model and 
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establishing its connection with the kdV hierarchy In the general case the inte- 
grable structure underlying the 0{n) model should likewise be encoded in the auxiliary 
function. The structure is well known when u is rational where the kdV„ hierarchies 
appear but the precise translation from matrix model variables to continuum time vari- 
ables is lacking. One might also ask whether there are integrable hierarchies associated 
with the continuum theories corresponding to non-rational values of z/, not to mention 
imaginary values of u. 
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